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TECHNICAL NOTE D-582 


PERFORMANCE CHARTS FOR MULTISTAGE ROCKET BOOSTERS 
By John S. MaeKay and Richard J. Weber 


SUMMARY 

Charts relating the stage propellant fractions are given for two- 
and three-stage rockets launching payloads into nominal low-altitude 
circular orbits about the earth. A simple method is described for ex- 
tending these data to higher orbit or escape missions. Various combina- 
tions of stages using RP - liquid- oxygen and hydrogen - liquid-oxygen 
propellants are considered. However, the results can be generalized with 
little error to any other propellant combination. 

Examples are given to illustrate how the charts may be applied to 
the preliminary design of booster vehicles. 


INTRODUCTION 

Optimal design of a multistage rocket vehicle for any particular 
mission involves a complicated compromise among such factors as mass- 
ratio requirements, structural weights, guidance, operational charactei - 
istics, and so forth. Of these, estimation of the needed stage mass 
ratios (or, equivalently, the stage propellant fractions) is ordinarily 
the primary design requirement, inasmuch as the propellant comprises the 
greatest volume and weight of most chemically fueled vehicles. 

Simplified methods for such estimations are usually' inadequate for 
vehicles operating in the vicinity of the earth because of the need to 
account properly for aerodynamic effects, gravity losses, varying thrust- 
weight ratios, and so forth. On the other hand, accurate calculations 
for each particular case are quite time-consuming because of the non- 
linearity of the differential equations of motion and the necessity for 
selecting an efficient flight path. To fill the need for a rapid, yet 
reasonably accurate technique, this report presents a group of charts 
for estimating the stage propellant fractions of booster vehicles that 
carry out orbital or escape missions. The accuracy of these generalized 
charts is adequate for preliminary design purposes. 

The two missions for which results are specifically given are flights 
under continuous power from the earth's surface into circular orbits at 
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altitudes of 75 and 150 nautical miles. As discussed later, optimum 
trajectories for more energetic missions usually pass through approxi- 
mately orbital conditions. Since the required superorbital man euvers 
can be closely approximated by impulsive- thrust equations, the presented 
data can be easily extended to cover any higher orbit or escape missions. 

The propellants considered are RP - liquid-oxygen and hydrogen - 
liquid-oxygen, used both throughout all stages and in various combinations. 
A simple method is given for modifying the data for propellants having 
other specific impulses. A range of thrust-weight ratios is covered that 
is typical of liquid- fueled rockets. 

SYMBOLS 

A area, sq ft 

a tan 

b d(tan \|r)/dt, sec -1 

p 

D drag coefficient, based on cross-sectional area 
D drag, lb 

F thrust, lb 

g local acceleration due to gravity, ft /sec*" 

g c veight-to-mass conversion factor, 32.174 lb/slug 

h altitude, ft 

I specific impulse, lb-sec/lb 

P ambient pressure, lb/sq ft 

R radius, ft 

r^ stage weight ratio, (Wq/w^.)^ 

t time measured from ignition of i^* 1 stage, sec 

V velocity, ft/sec 

AV characteristic velocity increment, ft/sec 

W weight, lb 
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w_ _■ propellant fraction, (W^/Wq ) j_ 

p, 1 r 

x x- coordinate of fig* 16, ft 

y y-coordinate of fig. 16, ft 

a angle of attack, deg (see fig. l) 

(3 angle between velocity vector and a normal to the radius vector 

(see fig. 1), rad (unless otherwise specified) 

Q angle between instantaneous and initial radius vector, (see fig. 

l), rad (unless otherwise specified) 

li gravitational constant, 1.4077X10 16 , ft 3 /sec 2 

p ambient density, slugs/cu ft 

angle "between thrust vector and x-axis, deg (see fig. 16) 

Subscripts: 

cs cross-sectional 

d evaluated at design altitude (h d = 20,000 ft) 

E evaluated at surface of earth 

e evaluated at rocket nozzle exit 

en engine 

f evaluated at turnout of i"^ stage 

i stage number 

L payload 

p propellant 

s structure 

t total 

0 evaluated at stage ignition (t = 0) 

Superscripts: 

— evaluated at t = Wq/w 
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relates to Impulsive velocity addition 
conditions after change of specific impulse 


ANALYSIS 

Assumptions 

tp 

The following assumptions were made in the calculations. Unless 
otherwise specified, the discussion relates to three- stage vehicles. 

Flight path . - As shown in figure 1, the vehicle takes off and 
climbs vertically for 10 seconds and is then turned for 10 seconds at a 
constant angle of attack a. A gravity turn (a = 0) is followed there- 
after until burnout of the first stage. During the next two stages the 
engine-thrust vector is directed according to the relation 


tan a - bt 

This linear- tangent ,r form of steering has been shown to result in approx- 
imately minimum propellant consumption for satellite launchings (refs. 1 
2, and 3). 

Drag coefficient . - A constant value of 0.4, which is typical of 
current vehicles, was used for the drag coefficient during first-stage 
flight. For the class of large boosters considered this corresponds to 
a value for ^D^cs/^0 1*6X10“^ square foot per pound and gives a con- 

servative estimate of the drag losses. No drag was Included during 
flight of the upper stages, which operate only at high altitudes where 
drag is negligible. (A typical first-stage burnout altitude for the 
cases considered is 180,000 ft.) 

Engine performance . - The propellant flow rate was maintained con- 
stant during each stage of flight. The specific impulse at sea level 
was taken as 267 seconds for RP-oxygen propellant and 358 seconds for 
hydrogen-oxygen; the assumed variation with altitude is given in figure 2. 
These values are typical of engines having chamber pressures on the order 
of 800 to 1000 pounds per square inch and nozzle expansion ratios of 15 
to 20. 


Equations of Motion 

The first-stage trajectories were calculated using a fourth-order 
Runge-Kutta numerical integration. Rotation of the earth was Ignored so 
that the data, as presented, apply to polar orbits. However, a method 
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for including earth rotation will be presented. The upper-stage linear- 
tangent paths were calculated with an approximate closed- form solution 
that is considerably faster than numerical integration. 

Details of the equations are given in appendix A. 


Procedure 


A number of first-stage trajectories were studied to obtain burnout 
conditions (altitude, velocity, path angle) as functions of v p^ p 
F/W 0 p, and a. These are shown in figures 3 and 4 for RP- oxygen and 
hydrogen- oxygen, respectively. Using values obtained from figures 3 and 
4 as initial conditions, the upper- stage trajectories were calculated 
for various combinations of F/w q in eacn stage and Wp in stages one 
and two. Automatic iteration was used on the steering parameters a 
and b and on third- stage burning time in order to achieve a circular 
orbit at 75 or 150 nautical miles. The output of each calculation is the 
required third-stage propellant fraction w p ^ 3. Each point was repeated 
for several values of a to find the minimum w p ^ 3. The calculations 
were performed on an IBM 704 digital computer and took approximately 25 
hours of machine time. 


Conditions along a representative flight path are shown in figure 5. 


RESULTS AN D DISCUSSION 
Propellant- Loading Charts 

The minimum values of w p ^ 3 have been plotted as functions of 
w^ -i and w-v> o. These basic data are given in figures 6 to 10 as indi- 

P, _L p, ^ 

cated in the following table: 


Figure 

Fuel 

in stage - 

Orbit altitude, 
naut. miles 

One 

Two 

Three 

e 

RP 

RP 

RP 

150 

7 

H 2 

H 2 

H 2 

150 

8 

RP 

RP 

H 2 

150 

Ci 

RP 

h 2 

h 2 

150 

10 

RP 
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For each case, all combinations of the following initial thrust-weight 
ratios are presented: 
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Stage 

F/W 0 

First 

1.20 

1. 30 

1.40 

Second 

1. 00 

1.25 

1. 50 

Third 

0. 50 

1.00 

1. 50 


Each point on the propellant-loading charts represents a three-stage 
vehicle capable of achieving the specified orbit. In general, higher 
orbits or escape energy is desired, and a correction must be made in the f 

propellant requirements. This is discussed in a later section. After 
making any necessary corrections, the determination of actual payload 
capabilities requires the incorporation of structural weights. Simpli- 
fied examples of such calculations are given in appendix B to illustrate 
how the charts may be applied. 

Thrust-weight ratio . - From the trajectory viewpoint high accelera- 
tion (thrust -weight ratio) is ordinarily thought to be desirable. This 
is because velocity is increased rapidly at low altitudes, and one saves 
the work required to lift propellants against gravity. However, the 
charts show that, for the 150-mile orbits, high acceleration tends to 
increase Wp^ 3 in some cases (e.g., fig. 6 (c)). The explanation lies 
in the requirement of the present analysis for continuous acceleration 
until a circular orbit is established. High-acceleration vehicles tend 
to build up too much velocity at low altitudes. Therefore, a less effi- 
cient boost path must be utilized in order not to exceed orbital velocity 
at burnout of the third stage. This problem is diminished if low orbits 
are desired, since the high velocities at low altitudes are then usually 
advantageous. Thus, it is seen that high thrust-weight ratios are more 
beneficial in the case of the 75-mile orbit. 

In any event, of course, no conclusions can be made about optimum 
thrust-weight ratios without including their effect on the weight of 
both the engines and the other thrust- sensitive structure. 

Drag coefficient . - A first-stage drag coefficient of 0.4 has been 
used in the charts. In order to indicate the importance of this param- 
eter, figure 11 shows the change in propellant consumption when the drag 
is taken as zero. It is seen that the results are quite insensitive to 
variations in the selected value. 


Extension of the Charts 

The charts as they are drawn apply to two specific missions, neither 
of which necessarily corresponds to the actual desired mission. This 
section describes how the charts may be applied to a much broader range 
of uses. 
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Except for very low orbit missions, flight paths incorporating a 
coast period are usually more efficient than those that continuously 
apply thrust. An example of this is shown in figure 12 for the 150-mile- 
orbit mission. The upper part of the figure reproduces the data from 
figure 8(c) for a w p 2 of 0.70. The lower part of the figure presents 
the result of optimum trajectory calculations that employed coasting 
between the second and third stages. It is seen that w p ^ 3 is improved 
by this procedure. (The crossover in the thrust-weight-ratio lines is 
also eliminated by including a coast period, for the reasons discussed 
in a previous section. ) The benefits of coasting are noticeable even 
for this comparatively low altitude of 150 nautical miles. They become 
more substantial as the final altitude is increased. 

Study of optimum trajectories that employ coast periods reveals the 
fortunate fact that the flight path prior to the start of coasting is 
usually nearly horizontal. Also, the velocity at that point is super- 
orbital. This means that very nearly optimum results would be obtained 
by achieving a circular orbit and continuing to accelerate until the 
proper velocity is reached, after which the vehicle coasts to its final 
altitude where it is injected into orbit. Calculations also show that 
the propellant consumption during both acceleration out of the low cir- 
cular orbit and injection into a higher energy orbit can be obtained 
with negligible error by the simple impulsive -thrust equations. The 
significance of these two facts is that, if data are given for achieving 
one initial low circular orbit, the total propellant consumption for 
achieving any higher orbit is readily derived. 

The choice of the best altitude for the initial circular orbit still 
remains. To show the effect of this selection, results are shown in 
figure 13(a) for a final circular orbit altitude of 500 nautical miles. 

In the range shown, the lower the initial altitude, the better. This is 
a result of the greater gravity loss when carrying fuel to high alti- 
tudes, the somewhat higher efficiency of departure from low orbits, and 
the interaction of altitude and thrust as previously discussed. 

Because of the presence of the earth's atmosphere, it is not prac- 
tical to lower the altitude indefinitely. Reducing the altitude requires 
the vehicle to achieve greater velocities while still in the denser por- 
tions of the atmosphere. The attendant greater drag penalties would 
eventually limit performance. Even more serious are the effects of aero- 
dynamic loads and heating. These structural factors lie outside the 
scope of this analysis, but it must be recognized that they exist and 
serve to limit the orbit altitude to some minimum value. For preliminary 
design purposes it is felt that the two orbit altitudes considered herein 
provide a range of trajectories which includes those of practical interest, 
(in addition, 150-mile orbits are high enough to be of some interest in 
themselves, such as for refueling or assembly purposes, yet are low 
enough to make continuous acceleration a reasonable practice if coasting 
should be undesirable, e.g., to avoid prolonged zero-g conditions or engine 

restarts . ) 
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The previous discussion has emphasized orbital missions, but the 
same argument applies to. escape missions. The only difference is that 
enough velocity is imparted when leaving the initial orbit to achieve 
escape energy. Of course, no injection into orbit at the end of coasting 
is required in this case. Figure 13(b) shows that the orbital— departure 
technique is a good approximation to the optimum direct-escape trajectory. 
(However, this good agreement is true only for the case of escape at 
approximately zero flight-path angle. Nonhorizontal escape may be occa- 
sionally desired, but specific trajectory calculations must then be made 
for that case. ) 

The charts discussed thus far have used three stages to gain the 
initial orbit. Another method of operation is to employ only two stages 
to reach the orbital condition. Two-stage data may be obtained from the 
charts (figs. 6 to 10) by extrapolating the curves down to the horizontal 
axis. For convenience, however, a chart specifically for this case has 
been prepared in figure 14 (which corresponds to an extrapolation of 
fig. 10). 

Procedure for extension . - As previously explained, the propellant 
consumption for more energetic missions is obtained by assuming the addi- 
tional velocity beyond the initial orbit is gained impulsively. It has 
been found that for top-stage thrust-weight ratios of 0.5 or higher, 
there is small error in this assumption. 

Figure 15 shows the propellant fraction w^ required for various 
characteristic velocity increments AV added beyond the no min al orbital 
condition. This is merely a plot of the simple impulsive relation 


AV = -Ig c ln(l - wp (1) 

Also indicated in the figure are the required values of total AV for 
various typical missions starting from initial circular orbits of 75 
and 150 nautical miles. Minimum- energy (Hohmann) transfer ellipses are 
assumed for the orbital missions. 

The values of w^ in figure 15 are based on the weight of the stage 
when in the initial orbit. If this stage (the third stage, say) has been 
already used in reaching the low orbit, its weight in the orbit is dif- 
ferent from its initial, fully loaded weight. The total stage propellant 
fraction after allowing for this difference is given by 

( v p,3)t = w p,3 + w p(! - W P , 3 ) (2a) 

where w p, 3 (which is the propellant used prior to the impulsive maneuver) 

is taken from the appropriate chart and (w p 3 ) t is based on the initial 
third- stage gross weight. ’ 
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Earth rotation . - Rotation of the earth was ignored when integrating 
the equations of motion. The presented data may therefore be interpreted 
as applying to launch azimuths of 0° or 180°. For nonpolar launchings, 
calculations showed that the charts may be applied with little error, if 
the additional AV used in figure 15 or equation (l) is reduced by the 
component of the earth's surface velocity lying in the launch direction. 
For example, an eastward firing from Cape Canaveral, Florida (latitude, 
28.5°) yields a benefit of 1350 feet per second. 

If no transfer to a higher orbit is to be made, the corrected w p ^ 3 
is given by 

( w p,3)t = w p,3 “ w p( X * w p,3) ( 2b ) 

where Wp corresponds to a AV equal to the surface velocity. 

Other propellants . - It was purely for the sake of convenience that 
the two propellant combinations considered herein have been designated 
RP - liquid oxygen and hydrogen - liquid oxygen. As far as the trajectory 
calculations are concerned, only the value of specific impulse is of 
direct interest. The charts may therefore be applied to any propellant 
combination yielding specific impulses identical to those of figure 2 . 

In addition, it :Ls possible to modify the presented data to obtain 
results for propellants having other values of specific impulse. This 
method depends on the fact that the stage AV is approximately independ- 
ent of the specific impulse of that stage. (This is not exactly true 
because varying the specific Impulse changes the variation of stage gross 
weight during flight. ) Assuming AVp does not change and applying 
equation (l) give the following expression for the new 1 after 
changes in the specific impulse of the stage (no change being made in 
the thrust-weight ratio): 

- £ - & - v I/r l w 

When using Ip in this equation, it should be an average value, which 
may be taken as one-third the sea- level value plus two- thirds the vacuum 
value. This gives 292 and 393 seconds for the RP and hydrogen data, 
respectively. 

Equation (3) is generally accurate to within 3 percent if Ip is 
varied, 2 percent if Ig is varied, and 1 percent if I 3 is varied. 

When several or all values of I are changed, the combined error seldom 
exceeds the error for variations of Ip. 
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Comparison with Mathematical Optimization 

Mathematical techniques for optimizing multistage vehicles have teen 
discussed ty many authors (e.g., refs. 4, 5, and 6). In general, these 
methods seek to find the stage propellant fractions that extremize a 
selected criterion (e.g., gross weight, empty weight, cost, etc.) subject 
to various constraints, such as turnout velocity, stage specific impulse, 
and so forth. 

For booster vehicles, the effects of drag, gravity, and turning 
losses are usually included ty either increasing the required turnout 
velocity or by depreciating the specific impulse of each stage in appro- 
priate amounts. Since, however, these losses are themselves functions of 
the propellant fractions, such corrections are mathematically invalid. 
Hence, it is not surprising that the numerical results of these techniques 
are usually in gross error. (A typical conclusion is that there should 
be no first stage, since the upper stages are more efficient. ) 

The deficiencies of the earlier techniques have teen remedied in a 
more recent paper (ref. 7). In this paper the trajectory and the vehicle 
are analyzed concurrently to find a true system optimum. However, the 
method requires first calculating the flight performance of a series of 
representative vehicles. Unless the optimum vehicle is fairly accurately 
known in advance, the amount of trajectory data that must be calculated 
is quite substantial - equivalent to generating one of the charts of the 
present report. 

If desired, the charts may be used in conjunction with the method of 
reference 7. However, direct calculations of the form shown in appendix 
B are probably simpler to carry out. 


CONCLUDING REMARKS 

Propellant- loading charts have been presented for multistage rockets 
launching payloads into nominal low-altitude circular orbits. Optimum 
trajectories for more energetic missions usually pass approximately 
through such orbital conditions. Since the required superorbital maneu- 
vers can be closely approximated by impulsive equations, the presented 
data for the nominal mission can be easily extended to cover any higher 
orbit or escape missions. Various combinations of stages using RP - 
liquid-oxygen and hydrogen - liquid-oxygen propellants are considered. 
However, the results can be generalized with little error to any other 
propellant combination. 

The data are based on trajectories that minimize propellant consump- 
tion. For real vehicles some of these optimum trajectories may engender 
excessive aerodynamic heating or structural loads. However, the 
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presented data are useful for preliminary design purposes even if the 
trajectories must be modified for some purpose. By incorporating struc 
tural weight assumptions, the charts permit the approximate optimum 
design of multistage chemical booster vehicles for either orbital or 
escape missions. The final optimization may be based on any desired 
criterion (e.g. , maximum payload, minimum cost, etc.); that is, the 
presented data do not restrict the optimization procedure. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, September 28, 1960 
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APPENDIX A 

METHOD OF CALCULATION 
Machine Trajectory Calculations 

First stage . - The first-stage trajectories were computed by numeri- 
cal integration of the equations of motion using the Runge-Kutta 
technique. 

In solving the equations of motion, the following simplifications 
were made: 

(a) Nonrotating earth 

(b) Constant drag coefficient (Cp = 0.4) 

(c) No lift. 


The two-dimensional equations in polar coordinates then are 


where 


. g c (F cos a - D) 


V = 


V3 = 


w 


V 


R E + h 


h » V sin 3 


\ _ V cos 3 
R E + h 


cos 3 + 


g sin 3 

g_F sin a 


W 


IW 


D = 


pV^CpA 


cs 


g = 


(Rg + h) £ 


The specif ic— impulse variations shown in figure 2 are given by 

1 - + ^ ( p a - p ) 


i OZJ 
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in which the design altitude is taken as 20,000 feet and no nozzle 
separation is assumed. 


During the integration process, the required values of P and p 
were found from a series of curve fits for the ICAO standard atmosphere. 

Each case was computed using a step size of 10 seconds and continued 

until the fuel consumption was such that w^ n » 0. 95. 

-tq x 

Second and third stages . - Neglecting drag and using tan \|r = a - bt, 
the equations of motion, written for the rectangular coordinate system of 
figure 16, are 


x + 


y + 


Hi = 

'Z 

IT 


g c F 


HI 

R 3 


W//(a - bt) 2 + 1 
gF(a - bt) 


= f. 


X 


W fj (a - bt )^ + 1 


- f. 




(Al) 


Assuming that the magnitude of R does not change during one stage of 
powered flight, equations (Al) become 


x + H2£ = f 
bS x 


y + ^ = f. 


(A2) 


R; 


o 


where Rq is the radius at the beginning of powered flight of the stage 
under consideration. 


By using the method of variation of parameters, equations (A2) can 
be solved in the form 


Y r = (y n + S) cos 03 t_p + I-ttt + B 1 sin cut 


(s 


(x 0 + U) cos ojt-o + l — - + E ) sin cut 


0) 


(A3) 
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where 



B 


E 



fy COS 03 1 dt 

— — 


f x cos crft dt 

~ — — 


> (A4) 


J 


Unfortunately, there Is no apparent way to integrate equations (A4) with- 
out some additional assumption. Con sider ing the magnitude of to, it^ can 

he seen that it will never exceed = 1. 2398X10" 3 radian per second. 

Consequently, if t < 200 seconds, which is a typical burning time for 
one stage, cut will always be less than 15°. For such small angles, it 
is reasonable to assume 


sin out = cut 
cos cut « 1.0 


Using equations (A5), equations (A4) can be integrated to give 


where 


S = V e [S(G - J sin \|r) + 

B = — (J sin If - G) 

U = cos ^ 

v e 

E = — J cos f 


ri + cos (tf - 1) cos \[r q" 
Ll + cos (\Jr q - \Jr) cos \|r . 


G = 


In 


[ 


(sin \|r^ 
(sin \|f 0 


+ l) cos >lr 0 ~| 
+ l) cos \j r J 


cos \|r^ - cos 
COS \[Tq cos 


(A5) 
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By using the 


the velocity 


fact that the method of variation 

• • • 

S cos cut + B sin cut = U cos cut + 

components can he written as 


of parameters assumes 
E sin cut = 0 


yf 


x 


f ~ 


-m(x 0 


+ S) sin cut^ 
+ U) sin cut^ 


+ cu 


+ cu 




+ B ) 

+ e) 


cos 


cos 


(X>tf 

cut^ 


Because of the assumptions made in equations (A5), the error of this 
approximate solution will increase with t^. An indication of how this 
error behaves is shown by the upper curve in figure 17, where the percent 
error in r^ is given as a function of tf 2 + ^f, 3 for several vehi- 
cles. The scatter* in the points shown is due to other differences in the 
vehicles chosen, but the essential feature to be noted is the rapid in- 
crease in error after tf 2 + tf^ 3 = 400 seconds (about 200 sec per 
stage). Consequently, it was necessary to use a modified solution for 
stages in which tf > 200. This modification consists of adding another 
term to the truncated series given in equation (A5). Thus, 


sin cut = cut 


(«*) 5 

3: 


cos cat = 1 - Mi! 

2: 


(A6) 


After substitution of equations (A6) into equations (A4), subsequent 
integration gives 


S 


V 

e 


tJ 




1 ] sin \Jr + 


(* * a •#(*¥- ") 


1 /co 2 t 2 


- 1 


of 


6b" 


2t 


2r \ . art tan \|r 
(tan t|t - 1) + 


2b 1 
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U = V. 


, T" 1 <X> 2 f 2 , . T 

t cos \|r \ — - — - 1 1 J + 


+ — (2 tan \|r - l) 
6b° 


E = V 


where 


(•¥->) 

! ]“ ' i ( ! • 1 * A - 


O _ 

, art tan t „ oo' 

2 — G 2 

2b 2b 


12b 


cos xgr 


CD 


co 2 ! 2 


) J .»h + i anj) G+ ^ 


cdH 

2b 2 


K = 


L = 


tan \|r f cos - tan \|tq cos 
COS \|Tq cos 

cos^q - cos 

’Z 

(cos 1 |Tq cos ty^) 


With these expressions for S, B, U, and E, the error is reduced to 
that shown by the lower curve in figure 17. 


Calculation of Impulsive Transfer Velocities 

Since energy and angular momentum are conserved along the path I- II 
(see fig. 18), 

A _ ji_ _ Ax _ _u_ 


Rn 


R 


II 


(A7) 

R I V I = R II V II ( A8 ) 

Using equations (A7) and (A8), Vj and Vjj can be solved for as 


V I " V c 


4 


2cp 


cp + 1 


II 


V 


-,^4 


tW + 1) 


where j - V^/ R I (25, 700 ft/sec for Rj - Rg = 75 naut. miles, and 
25,400 for Rj - Rg = 150 naut. miles) and cp = Rjj/Rj. Consequently, 
the required velocity increments (added at points I and II ) are 
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AV II ~ V C, II " V II \ I U ’ y V c P('P + i) J 

Note that as cp - AVj -> 0. 414 V c p which is the velocity increment 
required for escape* 


to 

i 

w 

o 
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APPENDIX B 


ILLUSTRATIVE EXAMPLES 

This section illustrates, by means of several highly simplified 
examples, how the charts may be applied to vehicle design. 

Example 1 . - Design a three-stage vehicle to place a 10, 000- pound 
payload into a 1000-nautical-mile orbit. Use the thrust-weight ratios 
and the simple weight relations given in the following table: 



Stage 


One 

Two 

Three 

Fuel 

HP 

h 2 

h 2 

F/W 0 

1.3 

1.0 

1.0 

W s /Wp 

.08 

.10 

.11 

Wen/F 

.010 

.012 

.012 


Assume that the vehicle is capable of withstanding the stresses incurred 
in trajectories that pass through orbital conditions at 75 nautical 
miles. As a start, pick w p ^ and w p ^ 2 as 0.7 and 0.6, respectively. 

From figure 10 (e), 3 = 0.395. From figure 15, the AV beyond orbit 

is 2850 feet per second. If we reduce this value by 1350 feet per second 
because of an eastward launch from Cape Canaveral, the corresponding 
propellant fraction from figure 15 (or eq. (l)) is 0.10. The total 
third- stage propellant fraction is given by equation (2) as 

w p ^ 3 = 0.395 + 0.10(1 - 0.395) = 0.456 
Now stun the component weights of the third stage: 


W 0,3 - W s 


W P + 


W en + W L 


' Vpf *K~) + w o(i^) ( t 2 ) 


+ W T 


Wt 


w„ 


Js\ 7 F \ / w en\ 
\ W p j - Vw 0 J ^ F j 


w, 


10^000 


°' 3 1 - 0.456(1.11) - (l.o) (0.012) 


= 20,750 lb 


Similarly, for the second and first stages (where the payload of each 
stage is the stage above it). 
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i 


20,750 


W 0, 2 1 - 0.6(1.10) - (1.0) (0.012) 


= 63,270 lb 


W 


63,270 


°^ 1 1 - 0.7(1.08) - (l. 3 ) (0. 010 ) 


= 273,900 lb 


The component weights, can now be computed also. 

To find the optimum staging, the previous calculations must be 
repeated for other v pjl and w p ^ 2 until W 0 ^ 1 (or any other selected 
parameter) is minimized. And finally, the entire optimization should be 
repeated for different values of thrust-weight ratios. 

Example 2. - Find the third-stage propellant fraction to escape from 
the earth if no coasting is permitted below 150 nautical miles. Use the 
data of example 1, but take the specific impulse of the hydrogen - 
liquid- oxygen engines as 400 seconds • 

From figure 9(ejj 3 = 0.430* Applying equation (3), the corrected 
values are 

w" , = 1 - (l - 0.430) 425 / 400 = 0.450 
v" = 1 - (1 - 0 . 60 ) 425 / 400 = 0.622 

p* 2 ' 

Figure 15 shows that an additional AV of 10,500 feet per second is 
required for escape. Using equation (l), 

10,500 = -(400) (32.174)ln(l - wp 
Wp = 0.557 

Finally, from equation (2a) 

w , = 0.450 + 0.557(1 - 0.450) = 0.756 
A weight analysis can now be made, as in example 1. 

Example 5 . - Determine the optimum method of staging for a 24-hour 
satellite mission. 

The mathematics for this problem is of the same nature as in the 
previous examples. The purpose of this paragraph is to point out the 
different methods of carrying out a mission that are possible using the 
charts. The following table lists how the various stages might be 
alternatively employed for this single mission: 
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Application 

Stage 

i 

ii 

in 

IV 

V 

VI 

VII 

VIII 

Boost into low orbit 
Depart from low orbit 
Inject Into high oribt 

1,2 

2 

2 

1,2 

2 

3 

CN3 

^ to to 

1 — 1 

i — i 

1, 2, 3 
3 
3 

1, 2, 3 

3 

4 

1, 2,5 
4 
4 

1, 2, 3 

4 

5 


Wo implication is intended that each of these eight alternatives is 
equally feasible. It is merely desired to emphasize the flexibility of 
the presented charts. 
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Figure 2. - Assumed variation of engine specific impulse with altitude. 
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(a) Velocity; F/W q 2 = 1.20. 

Figure 3. - F '.rat-stage burnout conditions. 1.^ = 292 (av.); c d a C3 Ao 1 = l.BxlO -4 . 
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(c) Velocity; F/W 


1.30. 


vexoci cy, r/w Q ^ 

Figure 3. - Continued. First-stage burnout conditions. = 292 ( av . ) ; Cj^ c , ; /Wq # \ ~ 1. 
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(b) Altitude and path angle; F/Wq ^ = 1.20. 

Figure 4. - Continued. First-stage burnout conditions. Ii = 393 (av.); CqA^/Wq ^ = 1.6*10" 
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(c) Velocity; F/W Q = 1.30. 

Figure 4. - Continued. First-stage burnout conditions. 1} - 393 (av.)j CdAcb/^O,! = 


1 .6X10 -4 . 
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(d) Altitude and path angle; F/W Q ^ = 1.30. 
Figure 4. - Continued. First-stage burnout conditions. 1^ = 393 (av. 
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(a) ¥/u 0fl = 1.20; F/W q ^ 3 = 0.50. 

Figure 6. - Three -stage propellant -loading chart. 150-Nautical 
mile orbit; ^ = 292 (av. ); I P = 305; I 3 = 305. 






First-stage propellant fraction, 

(b) F/W 0 ^ = 1.20} F/W^ 3 = 1.00. 

Figure 6. - Continued. Three-stage propellant-loading chart 
150-Nautica 1-mile orbit; I-, = 292 (av. ); I 2 = 305; I 3 = 30 
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First-stage propellant fraction, w p^i 

(c) F/W qa = 1.20; F/W 0 ^ 3 = 1.50. 

Figure 6. - Continued. Three-stage propellant-loading chart. 
150-Nautical-mile orbit; Ip = 292 (av. ); I 2 = 305; I 3 = 305. 
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Figure 6. - Continued. Three-stage propellant -loading chart 
150-Nautical-mile orbit; I-, = 292 (av. ); I P = 305; I 3 = 30 
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(g) F/W 0 ^ = 1.40; F/Wq^ = 0.50. 

Figure 6. - Continued. Three -stage propellant -loading chart 

150-Nautical-raile orbit; I-i = 292 (av. ); I 2 = 305; I 3 = 30 
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First-stage propellant fraction, w p^i 

(h) F/W 0 , = 1.40; F/W 0 ^ 3 = 1.00. 

Figure 6 . - Continued. Three-stage propellant -loading chart 
150-Nautical-mile orbit; I-. = 292 (av. ); In = 305; I 3 = 30 
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• 64 * 68 .72 . 76 


First-stage propellant fraction, Vp ^ 

(i) = 1.40; F/Wq^ 3 = 1.50. 

Figure 6. - Concluded. Three-stage propellant-loading chart. 
150-Nautical-mile orbit; = 292 (av. ); I 2 = 305; I 3 = 305. 
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.48 .52 .56 .60 .64 .68 .72 .76 

First-stage propellant fraction, p 


(a) F/W 0>1 = 1.20; F/w 0 3 = 0.50. 

Figure 7. - Three-stage propellant-loading chart. 150-Nautical-mile orbit; 1^ = 393 
(av. ); I 2 = 425; I, = 425. 
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.48 .52 .56 .60 .64 .68 .72 .76 

First-stage propellant fraction, 

(to) F/W 0>1 = 1.20; F/W 0?3 = 1.00. 

Figure 7. - Continued. Three -stage propellant -loading chart. 150 -Nautical -mile orbit; 
- 393 (av. ); Ig = 425; I 3 = 425. 
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First-stage propellant fraction, W P,1 


(c) F/W 0>1 = 1.20; F/W 0;3 = 1.50. 


Figure 7. - Continued. Three -stage propellant -loading chart. 

I-L = 393 (ay. )j I 2 = 425$ I3 - 425. 


150 -Nautical -mile orbit; 
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.48 .52 .56 .60 .64 .68 .72 .76 

First-stage propellant fraction, Wp^p 


(d) F/W Qj1 = 1.30; F/Wq 3 = 0.50. 

Figure 7. - Continued. Three-stage propellant-loading chart. 150-Nautical-mile orbit; 
I x = 393 (av. ); I 2 = 425; I 3 = 425. 
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Figure 7. - Continued. Three -stage propellant -loading chart. 150 -Nautical -mile orbit; 
I x = 393 (av. )$ I 2 « 425; I 3 - 425. 






(f) F/W 0>1 = 1.30} F/W 0>3 = 1.50. 


Figure 7. - Continued. Three-stage propellant-loading chart. 150-Nautical-mile orbit} 
!]_ = 393 (av. )} I 2 = 425} I 3 = 425. 
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(g) = 1-40; f/w 0>3 = 0.50. 

Figure 7. - Continued. Three-stage propellant-loading chart. 150-Nautical-mile orbit; 
I x = 393 ( av. ) ; I 2 = 425; I 3 = 425. 
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Third-stage propellant fraction, v p^3 
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(i) F/W 0jl = 1-40; F/W 0)3 = 1.50. 

Figure 7. - Concluled. Three-stage propellant-loading chart. 150-Nautical-mile orbit; 
I x = 393 (av. ); l ? = 425; I 3 = 425. 
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.60 .64 .68 .72 .76 .80 

First- stage propellant fraction, Wp ^ 

(a) F/W q x = 1.20; F/W q 3 = 0.50. 

Figure 8 . - Three -stage propellant-loading chart. 150-Nautical- 
mile orbit; = 292 (av. ); Ig = 305; I 3 = 425. 
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First-stage propellant fraction, w p^i 


(c) F/W 0 p = 1.20; F/W 0 ^3 = 1.50. 

Figure 8 . - Continued. Three-stage propellant-loading chart. 
150-Nautical-mile orbit; Ip = 292 (av. ); I 2 = 305; I 3 = 425. 
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(e) F/W qa = 1.30; F/W q ^ = 1.00. 

Figure 8. - Continued. Three-stage propellant -loading chart. 

150-Nautical-mile orbit; !•}_ = 292 (av. ); Ig = 305; I3 = 425. 
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First-stage propellant fraction, w n 

(g) F/W 0 ^ = 1.40; F/w q ^ 3 = 0.50. 

Figure 8. - Continued. Three-stage propellant -loading chart 
150 -Nautical -mile orbit; 1-^ = 292 (av. ); Ig = 305; 1^ = 42 
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(h) F/W 0 ^ = 1.40; F/W q ^3 = 1.00. 

Figure 8. - Continued. Three-stage propellant-loading chart 
150-Nautical-rnile orbit; I-, = 292 (av. ); I ? = 305: I, = 42 
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(i) F/Wq^ = 1. 40; F/Wq^ - 1.50. 

Figure 8. - Concluded. Three-stage propellant-loading chart. 
150-Nautical-mile orbit; = 292 (av. ); I 2 = 305; I 3 = 425. 
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(a) P/W 0 -l = 1.20; F/W q ^ = 0.50. 


Figure 9. - Three -stage propellant -loading chart* 150-Nautical- 
mile orbit} 1^ = 292 (av.)j !<? = 425} I 3 = 425. 
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First -stage propellant fraction, w p 

(b) F/W 0 ^ = 1.20; F/W 0 ^ 3 = 1.00. 

Figure 9. - Continued. Three-stage propellant-loading chart. 
150-Nautical-mile orbit; I ± = 292 (av. ); I 2 = 425; I 3 = 425. 
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(d) F/Wq^ = 1.30; F/W q ^ 3 = 0.50. 


Figure 9. - Continued. Three-stage propellant -loading chart. 
150 -Nautical -mile orbit; I 1 = 292 (av. ); I 2 = 425; I 3 = 425. 
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(h) F/Wq^-l = 1.40} J?/V Qf5 = 1.00. 


Figure 9. - Continued. Three -stage propellant -loading chart. 
150-Nautical-mile orbit; 1-^ = 292 (av. ); I 2 = 425; I 3 = 425. 
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.60 .64 .68 .72 .76 .80 

First-stage propellant fraction, w , 


(a) F/W q 1 = 1.20; F/W q 3 = 0.50. 

Figure 10. - Three-stage propellant -loading chart. 75-Nautical- 
mile orbit; I ± = 292 (av. ); I 2 = 425; I 3 = 425. 
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First-stage propellant fraction, w p^i 

(b) F/Wq^-l = 1.20; F/W 0#3 = 1.00. 

Figure 10. - Continued. Three -stage propellant -loading chart 

75-Nautical-mile orbit; = 292 (av. ); I2 = 425; I3 = 425. 
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•60 .64 .68 .72 .76 .80 

First-stage propellant fraction, w 

ir J 

(c) P/W 0 = 1.20} F/w Q 3 = 1.50. 


Figure 10. - Continued. Three -stage propellant -loading chart. 
7 5 -Naut i c al -mile orbit} = 292 (av. )} I 2 = 425} I 3 = 425. 
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First-stage propellant traction, w p^p 

(d) F/W q , = 1.30; F/w q ^ 3 = 0.50. 

Figure 10. - Continued. Three-stage propellant-loading chart 
75 -Nautical -mile orbit; Ii_ = 292 (av. ); I 2 = 425; I 3 = 425. 
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(e) F/Wq^ - 1.30$ F/Wq ^ 3 = 1.00. 


Figure 10. - Continued. Three-stage propellant -loading chart. 
75 -Nautical -mile orbit; 3^ = 292 (av. )$ I 2 = 425$ I 3 = 425. 
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First-stage propellant fraction, w p^i 

(f) F/W 0 ^ = 1.30; F/W 0 ^ 3 = 1.50. 

Figure 10. - Continued. Three-stage propellant-loading chart 

75-Nautical-mile orbit; In = 292 (av. ); I 2 = 4=25; I 3 = 425. 
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(g) F/W 0 ^ = 1.40} F/W q 3 = 0.50. 


Figure 10. - Continued. Three-stage propellant -loading chart. 
75-Nautical-mile orbit; I ± = 292 (av. ); I 2 = 425; I 3 = 425. 
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(h) F/W 0 ^ = 1.40} F/W 0;3 = 1.00. 

Figure 10. - Continued. Three-stage propellant -loading chart. 

75-Nautica] -mile orbit} Ip = 292 (av. )} I 2 = 425} I 3 = 425. 
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(i) F/Wq^ = 1.40; F/W Q , = 1.50. 


Figure 10. - Concluded. Three -stage propellant-loading chart. 
75-Nautical-mile orbit; 1 ^ = 292 (av. ); I 2 = 425; I 3 = 425. 
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First-stage propellant fraction, ± 

Figure 11. - Effect of drag on third-stage propellant fraction 
for various first-stage propellant fractions. w p^2 = 0.50; 

F/W 0 ! = 1.4; F/W 0 ^ 2 = 1”. 0; F/W 0 ^ 3 = 1.0; I ± = 292 (av. ); 

I 2 = 305; = 425; 150 -nautical -mile orbit. 
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•60 • 64 .68 .72 .76 .80 


First-stage propellant fraction, w 

Figure 12. - Comparison of continuous -‘burning ascent to 150- 
nauti cal -mile orbit with use of optimum coast period between 
second and third stages. F/w 0 p = 1.20; F/W q ^ = 1.50; 

W P,2 = °* 70 > Ip = 292 (av. ); I 2 = 305; I 3 = 425. 
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Initial orbit altitude, naut. miles 
(a) Transfer to 5 00 -nautical -mile circular orbit. 



(b) Escape missions. 


Figure 13. - Effects of impulsive thrust assumption and initial 

orbit altitude. F/Wq ^ = 1.40; F/Wq^ g = 1* ^A^O, 3 = -*-• 

w , = 0.60: v „ = o! 70; I, = 292 (av. ): Io = 305; I 3 = 425. 
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(a) F/W 0 x = 1.20. 


Figure 14. - Two-stage propellant -loading chart. 75-Nautical- 
mile orbit; 1^ = 292 (av. ); Ig = 425. 
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First -stage propellant fraction, 

(b) F/W 0 p = 1.30. 

Figure 14. - Continued. Two-stage propellant -loading chart. 
75-Nautical-mile orbit j Ip = 292 (av. )j Ig = 425. 
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First-stage propellant fraction, w p ^ 


(c) F/W 0 ^ = 1.40. 

Figure 14. - Concluded. Two-stage propellant -loading chart. 
75 -Nautical -mile orbit; I-^ = 292 (av. ); I 2 = 425. 
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Figure 17. - Variation of third-stage weight-ratio error with total upper 
satellite boosting vehicles. 150 -Nautical -mile orbit. 





